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ON THE STABILITY OF CYCLES BY DELAYED 
FEEDBACK CONTROL 
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Dedicated to Professor Yuriy Kryakin on the occasion of his 60 th birthday 


Abstract. We present a delayed feedback control (DFC) mecha¬ 
nism for stabilizing cycles of one dimensional discrete time systems. 
In particular, we consider a delayed feedback control for stabilizing 
T-cycles of a differentiable function / : R —>• R of the form 

x(k + 1) = f(x{k)) + u(k) 

where 

u(k) = (a 1 -l)f(x(k))+a 2 f(x(k-T))-\ - \-a N f(x(k-{N-l)T)) , 

with a± + ■ ■ ■ + a^v = 1. Following an approach of Morgiil, we 
construct a map F : R T+1 —> K T+1 whose fixed points correspond 
to T-cycles of /. We then analyze the local stability of the above 
DFC mechanism by evaluating the stability of the corresponding 
equilibrum points of F. We associate to each periodic orbit of / 
an explicit polynomial whose Schur stability corresponds to the 
stability of the DFC on that orbit. An example indicating the 
efficacy of this method is provided. 


1. Introduction 

The control of chaotic systems is a subject of considerable interest in 
a number of disciplines, in particular engineering, physics, and mathe¬ 
matics. In the fundamental paper [I], Ott, Grebogi, and Yorke observed 
that chaotic systems frequently contain unstable periodic orbits that 
may be stabilized by small time-dependent perturbations. Subsequent 
papers such as have explored specific control mechanisms for sta¬ 
bilizing chaotic systems. One method of control of particular interest 
is the delayed feedback control (DFC) scheme introduced by Pyragas 
in pfj. The control in the Pyragas scheme is essentially a multiple of 
the difference between the current and one period delayed states of the 
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system. Primary advantages of this scheme include the facts that the 
control term vanishes if the system is already in a periodic orbit and 
that the control term tends to zero as trajectories approach a given 
periodic orbit. The DFC control mechanism enjoys many applications 
ranging from the stabilization of the modulation index of lasers to the 
suppression of pathological brain rhythms mm- 

In spite of its relative simplicity and broad range of application, the 
stability analysis of the DFC mechanism remains a difficult issue. A 
paper of particular interest to us in this regard is one of Morgiil. In 
[7], Morgiil considers the one-dimensional discrete time system 

x{k + 1) = f{x(k)) + u(k) , 

with k G Z being the time index and / : R. —> R an appropriately 
differentiable function. We suppose that / has a (possibly unstable) T- 
periodic orbit T t = {x* 0 ,x* T _ x }, where f(x* mod T ) = f{x* +1 mod T ). 
A DFC control used to stabilize the orbit Tt is 

u(k) — K(x(k) — x{k — T)). 

As shown by Morgiil, we may analyze the local stability of this control 
by considering the auxiliary function G : R r+1 —y R T+1 defined by 
G(zu • • •, z T +0 = (z 2 , ■ • •, Z T+ 1 , f(zr+ 1 ) + K{zt+i - 2 i)). We define 
F : R T+1 —>■ R T+1 by F — G T (the composition of G with itself T 
times.) Observe that T' t := {xq, ..., Xy_j, x^} is a fixed point of F. 
The stability of Tt under the control mechanism is equivalent to the 
stability at T' t of the system x{k + 1) = F(x{k)). The latter may be 
analyzed by computing the Jacobian of F at T' t . The Jacobian of F has 
a characteristic polynomial p( A), and the T-cycle Tt is exponentially 
stable under the control provided that p( A) is Schur stable, i.e. all 
of its eigenvalues lie inside the unit disc of the complex plane. In [7] 
Morgiil explicitly provided the calculation of the above characteristic 
polynomial. Setting a\ = /'(xq), ... ,cit = /'(^r-i); the characteristic 
polynomial of the Jacobian is of the form 

p{ A) = + ct^ t + • • • + ci A + Co 

where c 0 = —(—1 ) T K T , ct = — riiLi( a * + AT), and 

ct-«=-(-i ) i K< t £ - t n (<x+^)- 

*1=1 *2=»i+l U=it- l+l i=1 

ig 
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Motivated by this previous work, we wish to consider now a control 
that takes into account a deeper prehistory of the output values of a 
function. In particular, we wish to consider a control of the form 

u{k) = (cii - 1 )f{x(k)) + a 2 f(x(k — T)) -\ -f a N f(x(k - (N - 1)T)) , 

where ai + • • • + ajv = 1, that takes into account not only the value of / 
at Xk-r but also Xk- 2 r, Xk- 3 T, ■ ■ ■, Xk~(N~i)T- The reason for consider¬ 
ing a control of this form is that, rather than having only one parameter 
K that may be modified in our attempt to provide stability, we have 
a set of parameters cq,..., a at at our disposal that may be adjusted to 
provide a more robust control mechanism. 

Proceeding analogously to the ideas of Morgiil, we associate to the 
above control a map G : R T ( JV ~ 1 )+ 1 —y ]gT(v-i)+i defined by 

G(x i,..., xt(at-i)+i) = 

(x 2 ,x 3 ,..., xt(n- i)+ij «i/(^T(iv-i)+i)+« 2 /(*Pr(iv— 2 )+i) T■ • -+aNf{x i)). 

We dehne F = G T , and the stability of the above control on the 
cycle T,t may be ascertained by the location of the roots of the charac¬ 
teristic polynomial of the Jacobian of F at (xq,xI, ... ,x^ N _^) where 
for convenience we set x* = x* mod T . (We will see that the roots of 
this characteristic polynomial are independent of which point in the 
T-cycle xq corresponds, a fact that is not a priori obvious.) 

This approach is all well and good but is of little value if one cannot 
explicitly compute the characteristic polynomial of the above Jacobian. 
The prospect of doing may strike one as rather intimidating, especially 
seeing the result of Morgiil’s calculation indicated above when one is 
considering a control with only one free parameter. The primary pur¬ 
pose of this note is to show that the desired characteristic polynomial 
p( A) may not only be computed but moreover has surprisingly simple 
form 

PO) — A (W_1 ’ r+1 - , 

where /i = /'(xq) ■ ■ ■ /'(^t-i) an d 

q'(A) — aiA T ■ ■ ■ T &N —iA T cijv • 

The derivation of the polynomial p( A) was, at least to us, decidedly 
unobvious, and it seems appropriate for it to be available in the litera¬ 
ture. This is the subject of the next section. In the subsequent section 
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we will indicate an application of this new control and advantages it 
lias over a control only taking into account a prehistory of length T. 

2. Derivation of the polynomial p( A) associated to the 

CONTROL 

We begin by providing the Jacobian associated to the control 

u(k) = (ai - l)f{x(k)) + a 2 f{x(k — T)) H-f a N f(x(k — (N — 1 )T)) 

at a point x* = (xj$, x*,..., x^ n _^) G R t ^ jv_1 ) +1 , where {xq, ..., x^_ x } 
forms a T-cycle of / and x * +1 = /(x*). For convenience we assume / 
is everywhere continuously differentiable and set jij = /'(x*) and 
/i = fiifi 2 ■ ■ ■ Ht- Recalling that G(x i,... ,xt(n- i)+i) = 

(x 2 , X 3 , . . . , Xt(AT- 1 )+ 1 , dlf{xT(N-l)+l) + ^2f(xT(N-2)+l) + - ' , 

we see by means of induction and application of the chain rule that the 
Jacobian of F — G T is given by the (N — 1 )T + 1 by (N — 1 )T + 1 
matrix J whose entries are given by 

'l if j = i + T 

i—{N—2)T—l 

■ \ _ a v ( ' N 2T ^ Hk if j = s mod T, i > (N — 2 )T + 2, and 

k=s 

(j — 1) mod T < {i — {(N — 2 )T + 2) 
^0 otherwise . 


The above expression is very nonintuitive, and as a help to the reader 
we provide an example of such a matrix when N = 3 and T = 3. 


( 0 

0 

0 

1 

0 

0 

0 \ 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

03/G 

0 

0 

02/U1 

0 

0 

ai/ii 

aid3pi(j,2 

a 3 /i2 

0 

O1O2/41/0.2 

02P2 

0 

a\/iiii2 

\ a\a3HiH2H3 

0103/^2/43 

03/13 

afa2/ii/i2/43 

01O2/O2/43 

02/43 

d\n 1/12/43 / 


Figure 1. Jacobian Matrix for N— 3, T —3 


The following lemma and subsequent theorem constitute the heart 
of the mathematical content of this paper. 
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Lemma 1. Let J be the (N — 1 )T + 1 by (N — 1 )T + 1 matrix given 
above. The eigenvalues of J are precisely the roots of the polynomial 

p{ A) = A (7V ~ 1)t+1 - Oi • ■ •/i T )(aiA Ar_1 H-f a N ~iX + a N ) T . 

Proof. We begin by quickly dispatching the T — 1 case. Observe that 
when T = 1 the Jacobian matrix takes on the form 

/ 0 1 0 0 ••• 0 \ 

0 0 1 0 ••• 0 

0 0 0 1 
\ GpvP-l GpV—I/O a N-2pl ' ' ' a lh 1 / 

This is a classical “companion matrix” (see, e.g. 0 ) and its eigenvalues 
are the roots of the polynomial 

p{ A) = \ N - pi(a\\ N ~ l + a 2 X N ~ 2 H-h ajv) , 


as desired. 

In the following we assume that T > 2. Let A be an eigenvalue of J. 
Then J(x i,..., X(at_i)t+i) T = A(xi, ..., X(at-i)t+i) T for some nonzero 
vector (xi,..., X(tv_i)t+i) t . (Here to avoid confusion we are using the 
upright “T” to denote the transpose of a matrix.) For our convenience, 
we set 


q(X) — u_/v T otv— iX + • • • + aiA^ 1 


and 


r(A) = q( A) — ai X N 1 . 

As J(efc) = ek+r for k — 1,..., T — k, we immediately have that A must 
satisfy the following simultaneous equations for nonzero (x±,. .. , Xt) T - 


' xiHiq(X) = X n_ 1 x 2 
Xiai/ii/i 2 g(A) +x 2 /ji 2 r( A) = A 7 V_ 1 x 3 
< : 

xiaf~ 2 /ii • • • fiT-iq(X) + (x 2 ai _3 /i 2 • • • /ir-i H-h x T -iPT-i)r(X) = X N ~ l x T 

''' TrqiX) + (x 2 ai~ 2 ii 2 ■ ■ ■ p T H-b x T piT)r( A) = A^xi . 
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This holds hold for nonzero (aq,..., Xt) t if and only if the determinant 
of the T by T matrix 

° \ 
0 

0 

af~ 2 /ii •'' /ir-ig(A) af _ 3 /r 2 • • • /r T -ir(A) • • • p T -ir(X) -A ^ -1 

V af ~Vi' • • A*r?(A) - X N ol~ 2 ii 2 Hz ■ ■ ■ Hrr(X) ■ ■ ■ ai/i T _i/i T r(A) /i T r(A) / 


/ /^i9(A) 

ai/iiyU 2 g(A) 
a 2 /ii/i 2 /i3q , (A) 


—A 

/i 2 r(A) 

cq/^ 2 /w(A) 


0 0 

-A^' 1 0 

/i 3 r(A) -A^ 1 


equals 0. It suffices to show that this determinant is equal to the 
polynomial —p( A) = — \( N ~ 1 ) T + 1 + (/q • • • /i T )(g(A)) r . 

We proceed by induction on T. Note that 


det 


hM A) -A^ -1 \ 

ai/ii / n 2 g(A) - \ N /i 2 r(A) y 


= ^ 2 (g(A) 2 - A 2 ^ 1 . 


So the T = 2 case holds. 

Suppose now that the T — (k — 1) case holds, i.e. 


/ ^ffi(A) 

ai^iAt 2 g(A) 


det 


-A^" 1 0 

/U 2 r( A) -A^ -1 


V 4 Vi • • • /ifc-ig(A) - X N a\ 2 /i 2 • •-/ifc-ir(A) 


\ 


-A^- 1 
Hk-ir(X) ) 


= H1H2 • • ■// fc _i(g(A )) fc - 1 - A (iV_1)(fe_1 ) +1 

Then, taking advantage of the fact that 


/ Xi 


det 


0 \ 


x 2 


5 
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we have that 

/ hi<?(A) 

aihih2<?(A) 


-A^" 1 0 

/12F(A) -A^ -1 


det 


\ 


-A^- 1 

\ a \~Vi ■ • • Hk-iq(X) a\~ 2 H2-■ ■ Hk-ir{X) ••• Hk-ir(X) / 


= hi/i 2 • • •/i fc -i(g(A )) fc - 1 - a^-W *- 1 ^ 1 + A 7 V (A Ar - 1 ) fc - 2 
= Hi • • •/ifc_i(g(A )) fc_1 . 


So 


/ hi^(A) 

ai/ii/i 2 g(A) 


det 


-A^- 1 

/i 2 F(A) 


0 

-A^- 1 


IV-1 


\ a\n 1 • • • /ifcg(A) - A w a* V 2 ■ ■ • /U fc r(A) 


-A 

^fcr(A) / 


/ m(A) 

aiA t iA t 29(A) 


= /ifcF(A) det 


-A^- 1 

/i 2 r(A) 


0 

-A^" 1 


\ 


V a, Vi • • • /Xfe-ig(A) a, >2 • • • hfc-i r (A) 


-A^" 1 
/i fc _ir(A) / 


/ /d<?(A) 

ai/ii/i29(A) 


S-A^ 1 det 


-A^" 1 

/i 2 r(A) 


0 

-A^" 1 


a? Vi • ■ • Hk-2q(X) a\ 3 H2-■ ■ Hk-2r(\) 
\ a*Hi ■ ■ ■ hfc<?(A) - \ N a* - V 2 • • • /i fc r(A) 

= /w(A) [hi • • •hfc-i(<?(A)) fc_1 ] 
+A JV ~ 1 


-A^- 1 
dlHk-lHk r {X) ) 

(induction) 
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= /ii ii'-WUiA" 1 (g(A)) 1- ' 1 (q , (A)) s '-A 1 -' 1 1)l!+1 

= / Jl ---^(9(A))‘-A( N - 1 > t+1 . 

Hence by induction we see the result holds. □ 

Theorem 1. Let J be the (N — 1 )T + 1 by (N — 1)T + 1 matrix given 
above. The characteristic polynomial of J is the polynomial 

p{ A) = A (Ar_1 ) T+1 - (hi ■ ■ ■ p. T )(ai\ N ~ 1 H-f a N -iA + a N ) T . 

Proof. The characteristic polynomial of J is of course given by det (A I — J), 
where I is the (N — 1)T + 1 by (N — 1)T + 1 identity matrix. Observe 
that this is a polynomial of degree (iV — 1)T +1 with leading coefficient 
1. Hence, p( A) will be the characteristic polynomimal of J, provided 
all the roots of p( A) are distinct, as by the lemma p( A) is a polynomial 
of degree (N — 1 )T + 1 with leading coefficient 1 and sharing the same 
roots as the characteristic polynomial of J. 

We now deal with the case that p( A) has a root of higher multiplicity. 

The idea we employ is that there should exist a minor perturbation of 
ai,... ,a,N that leads to a new polynomial with distinct roots, and the 
continuity of the characteristic polynomial of J in terms of a±,... ,On 
should then imply that p( A) is indeed the desired characteristic poly¬ 
nomial of J. 

Proceeding formally along these lines, we fix ai,..., ajv, and for z e 
C we define the matrix J zw as the (N — 1)T +1 by (N — 1)T +1 matrix 
whose entries are given by 

1 if j — % -\-T 

i—(N—2)T—l 

a z w jv_ i ini i <-ir 2T) - (s+i) n j - s m ° d t * - ( n - 2 ) t + 2 ’ and 

i>3) = ' ' T k= s 

(j - 1) mod T < (i - ((N - 2)T + 2) 

0 otherwise . 

where a z . w ,j = zaj for j = 1,..., N — 1 and a ZtWy N = a at + w. This 
is simply the previous Jacobian matrix J but where each a* has been 
replaced by za,i for i — 1,..., N — 1 and by + u>- Now, it suffices 
to show that there exists values of z arbitrarily close to 1 and w close 
to 0 for which J ZtW has no repeated eigenvalues. Equivalently, it suffices 
to show there exist values of z arbitrarily close to 1 and w close to 0 
for which the polynomial 

p z ,w{ A) = A (Ar_1)T+1 - (/ii • • • p T )(zai\ N l H-f za N - X \ + a N + w) T . 

has no repeated roots. To show this, we take advantage of the fact 
that the polynomial p ZjW ( A) will have a repeated root if and only if its 
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resultant R(p ZiW ,p' zw ) is zero. (See, e.g., [9].) Recall that the resultant 
R(f,g ) of the polynomials f{x) = a n x n + a n _ \X n ~ l + • • • + a^x + a 0 
and g{x) = b m x m + b m _ ix m_1 + • —h b\X + b 0 is the determinant of the 
Sylvester matrix 


( 


Q"n— 1 
dn 


^n —1 


a 0 


^n —1 


ao 


ao 


\ 


bm 


bm —1 
bm 


bm —1 
bm 


bo 


bm —1 


&0 


®n—1 


&0 


ao 


V 


• b m —i 


&0 / 


We then observe that the resultant R(p z , w ,p z w ) may be viewed as a 
holomorphic function in C 2 . Now, if p z>w ( A) had repeated roots for all 
z sufficiently close to 1 and all w close to 0, then R{p z , w p' zw ) would be 
identically zero for all (z,w) G C 2 . (See, e.g., [lOj H].) Setting z = 0 
and w = 1 — Otv, we would then have that the polynomial A( 7V-1 ) T+1 — 1 
had repeated roots, a clear contradiction. 

□ 


3. Stabilization of 1-cycles 

We now illustrate the effectiveness of the control 

u{k) = (ai - 1 )f(x(k)) + a 2 f(x(k -T))-\ -f a N f(x(k - {N - 1 )T)) 

in the case that T — 1 by showing that for any value g := pi < 1 the 
characteristic polynomimal of the associated Jacobian, namely 

p{ A) = X N — n(ai\ N 1 + • • • + a^v-iA + a^v) , 

has roots only in the unit disc for appropriately chosen N and ai,..., a^. 
We remark that the control will automatically be unstable for p > 1 
for any choice of N and a* with a\ + • • • + = 1 as the associated 

characteristic polynomial p{ A) would satisfy p( 1) = 1 — p < 0, thus 
failing the Jury test. Moreover, if 0 < p < 1 we may simply set 
ai — 0,... ajv-i = 0, oat = 1 for any N > 1 to achieve a polynomial all 
of whose roots lie in the unit disk. So it suffices to consider only the 
case that p < 0. 
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Let N be a positive integer. We set aq,..., on to all be 1/N. A root 
A of the above polynomial would then need to satisfy the equation 

° = A 7 V --^(A JV - 1 + "- + A + 1 ) . 

Now, if n — 0 the only solutions of the above equation is A = 0 
which of course lies in the unit disc. Moreover, the complex roots 
of a polynomial are continuous in the coefficients of the polynomial 
(see, e.g., p( 2 ]) and accordingly the solutions of the above equation will 
lie inside the unit disc for all values of p down to the greatest (but 
negative) value of 7 such that 

0 = A 7 V -^(A 7V - 1 + "- + A + l) 

has a root A that lies on the boundary of the unit disc. 7 is given by 
the equation 

— = — inf{9?(A _iV H— • + A” 1 ) : A G <90 and A(A _JV + • • • + A -1 ) = 0} . 
7 N 

Expressing values on the boundary of the unit disc as e ie for 9 G M, we 
have 


- = — inf{9f?(e- i7ve + - • • + e ~ ie ) :0GK and ^(e~ iNe + - • -+ e - i0 ) 
If N is odd, we may reexpress e~ lNe + • • • + e~ ie as 


= 0 }. 


-iNG 


H-he = e 


-iG —i (JV A 1) 9 f -i^-=^G 


H-+ e l ~ 


= e 


-iiEpie I 

sin ^ 


taking advantage of the closed form for the Dirichlet kernel. We see 
that the imaginary part of the above expression is zero only if 6 is 
either of the form for some integer k (for which sin(^) = 0) or of 
the form for some integer k (for which the imaginary part of the 
exponential term vanishes.) Exploiting the symmetry of the summands 
about the origin, we have that e~ lNe + • • • + e~ le = 0 if 9 = ^ and 
also that e~ lNe + • • • + e~ l9 = — 1 if 9 = . Hence 

-1 = 1 inf{9fj( e - i7ve +- • • +e ~ ie ) : 9 G M and $s{e~ m +- ■ -+e- ie ) = 0} . 

Accordingly, we see that by setting aj = 1/N for j = 1, ...,AT, the 
above control will be stable for all values of p in (-N, 1). In particular, 
regardless of how large the magnitude of /i, we see that by choosing N 
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sufficiently large the above control will be stable about the equilibrium 
point associated to fi. 

4. Conclusions and Future Directions 
We have seen that that the stability of the DFC control mechanism 
x(k + 1) = f(x(k )) + u(k), 

where 

u(k) = (ai - 1 )f(x(k)) + a 2 f(x(k — T)) -\ -f a N f(x(k - (N - 1)T)) 

with ai + • —h = 1, about a T-cycle may be ascertained by consid¬ 

ering the Schur stability of a polynomial of the form 

p( A) = A (7V-1)r+1 - ^(aiA^ -1 H-f a N -iX + a N ) T . 

As an application, we showed in the case of 1-cycles that this polyno¬ 
mial is Schur stable for appropriately chosen N and ai,..., a at for any 
given value of p, < 1. Moreover, in the 1-cycle case we have seen that 
this may be done where N is of the same order of magnitude as |/i|. 

A natural subsequent problem is the following: 

Problem: For any given integer T > 2 and p < 1, do there exist N 
and associated a\,... ,cin where ai + • • • + a at = 1 such that all of the 
roots of the polynomial p( A) lie in the unit disc? If so, what would be 
the smallest value of N for which this could be done? 

The reader may at this point be speculating whether or not in the 
T > 2 case one may simply choose uniform coefficents a\ = ■ ■ ■ = = 1/N 

and then try to prove the roots of p( A) all lie in the unit disc for suf¬ 
ficiently large N. We have analytic evidence that this cannot be done. 
However, Dmitrishin and Khamitova proved in [13] that the T = 1 case 
may be done with an improved and optimal value of N on the order of 
magnitude of p | 1 2 , where the associated a,j are given by the formula 

aj = 2 tan ———- • (1- 3 —\ ■ sin —j = 1,..., N . 

0 2(JV + 1) V A^ + iy N+l J ’ ’ 

Moreover, in the paper [ITj, the authors indicate that the T — 2 case 
may be done with a value of N on the order of magnitude of |/ii/U 2 |, 
and moreover that the order of magnitude in terms of the exponent 
of I/X 1 /X 2 1 is sharp. A natural conjecture at this point is that T-cycles 
associated to a multiplier Hi ■ ■ ■ pr may be stabilized by the above con¬ 
trol for appropriate values of ai, ..., ajy where N is on the order of 
magnitude of \p 1 ■ ■ ■ Pt\ t ^ 2 ■ The associated proofs of the above results 
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involve delicate estimates associated to trigonometric sums and more¬ 
over suggest that refined techniques in harmonic analysis may provide 
an avenue for affirmatively answering this problem, hence establishing 
a stability result for the above control regardless of the the length of 
the cycle T and the magnitude of the associated multiplier fi. This is 
a subject of ongoing research. 
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